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Abstract 

We show that the structural properties and phase behavior of a self-avoiding polymer chain on 
adhesive substrate, subject to pulling at the chain end, can be obtained by means of a Grand 
Canonical Ensemble (GCE) approach. We derive analytical expressions for the mean length of 
the basic structural units of adsorbed polymer, such as loops and tails, in terms of the adhesive 
potential and detachment force, and determine values of the universal exponents which govern their 
probability distributions. Most notably, the hitherto controversial value of the critical adsorption 
exponent <p is found to depend essentially on the interaction between different loops. The chain 
detachment transition turns out to be of the first order, albeit dichotomic, i.e., no coexistence of 
different phase states exists. These novel theoretical predictions and the suggested phase diagram 
of the adsorption-desorption transformation under external pulling force are verified by means of 
extensive Monte Carlo simulations. 
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Introduction - The manipulation of single polymer chains has turned recently into an im- 
portant method for understanding their mechanical properties and characterization of the 
inter.no.ecu.ar interaction, fl Q. Such m an,pu.ation is mainly triggered by the progress 
in atomic force microscopy (AFM) [3| as well as by the development of optical/magnetic 
tweezers technique |4|. This rapid development has been followed by theoretical considera- 
tions, based on the mean - field approximation [5], which provide important insight into the 
mechanism of polymer detachment from adhesive surfaces under the external force pulling. 
A comprehensive study by Skvortsov et al. |6| examines the case of a Gaussian polymer 
chain. We also note here the close analogy between the force detachment of adsorbed chain 
and the unzipping of a double - stranded DNA. Recently, DNA denaturation and unzip- 
ping have been treated by Kafri et al. 71 using the Grand Canonical Ensemble (GCE) 

n n 

approach |8| as well as Duplantier's analysis of polymer networks of arbitrary topology |9J. 
An important result for the properties of adsorbed macromolecule under pulling turns to be 
the observation [7j that the universal exponents (which govern polymer statistics) undergo 
renormalization due to excluded volume effects, leading thus to a change of the order of 
DNA melting transition from second to first order. In this Letter we use similar methods 
to describe the structure and detachment of a single chain from a sticky substrate when the 
chain end is pulled by external force. 

Single chain adsorption - Starting with the conventional (i.e., force-free) adsorption, we 
recall that an adsorbed chain is build up from loops, trains, and a free tail. One can treat 
statistically these basic structural units by means of the GCE approach where the lengths 
of the buildings blocks are not fixed but may rather fluctuate. The GCE-partition function 
is then given by 

A jv _ Vo(z) Q(z) 
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where z is the fugacity and U(z), V(z), and Q(z) denote the GCE partition functions of 
loops, trains and tails, respectively. The building block adjacent to the tethered chain 
end is allowed for by Vo(z) = 1 + V(z). The partition function of the loops is defined as 
U{z) = Yl™=i (A i 3- 2 ) n / nQ 5 where /X3 is the 3d connective constant and a is the exponent 
which governs surface loops statistics. It is well known that for an isolated loop a = 
1 — 711 « 1.39 [10|. We will argue below that a changes value due to the excluded volume 



interaction between a loop and the rest of the chain. The train GCE-partition function 



reads V(z) = Y^mLi (fJ>3 , wz) n /n 1 ~ ld = 2 whereby one assumes that each adsorbed segment 
gains an additional statistical weight w = exp(e) with the dimensionless adsorption energy 
e = e/ksT. Eventually, the GCE partition function for the chain tail is defined by Q(z) = 
1 + YunLi (l I 3 z ) n / n ^- F° r an isolated tail (3 — 1 — 71 ~ 0.32 [h]] but again the excluded 
volume interactions of a tail with the rest of the chain increase the value of f3. Using the 



generating function method [llj, S/v is obtained as = (z*)~ N where the pole z* is given 
by the condition V(z*)U(z*) = 1 so that the free energy is F = JzbTN In z* and the fraction 
of adsorbed monomers n = — dlnz* /din w. In terms of the so called polylog function, which 



is defined as &(a,z) = YlnLi z n /n a [12J], the equation for z* reads 

§(a,fi 3 z*)$(\,[i2wz*) = 1. (2) 

A nontrivial solution for z* in terms of w (or the adsorption energy e) appears at the critical 
adsorption point (CAP) w = w c where w c is determined from C(a)$(l —7^=2, [i2W c / ^3)) = 1 
as z* = 1/fi-s and ((a) is the Riemann function. In the vicinity of the CAP the solution 
attains the form 

z*(w) w [1 - A (w - wj 1 ^- 1 ^ 1 (3) 

where A is a constant. Then the average fraction of adsorbed monomers is n oc (e — 
ec) 1 '^ -1 ) -1 . A comparison with the well known scalin g re lationship n oc (e — e,,) 1 ^^ 1 where 
(j> is the so called adsorption (or, crossover) exponent [lOj suggests that 

= a - 1 (4) 

This is a result of principal importance. It shows that the crossover exponent 0, describing 
polymer adsorption at criticality, is determined by the exponent a which governs polymer 
loop statistics! If loops are treated as isolated objects, then a — 1 — 711 ~ 1.39 so that 
= 0.39. In contrast, excluded volume interactions between a loop and the rest of the chain 
lead to an increase of a and 0, as we show below. 

Probability distributions of loops and tails - How does the length distribution of polymer 
loops and tails close to the CAP look like? From the expression for U(z), given above, 
and eq.(ED we have Pi oop w (faz*) 1 w exp[-ci(e - e c ) 1 ^]/l 1+<1 '. This is valid only for 
e > e c since a solution for eq.© for subcritical values of the adhesive potential e does not 
exist. Nontheless, even in the subcritical region, e < e c , there are still monomers which 



occasionally touch the substrate, creating thus single loops at the expense of the tail length. 
The partition function of such a loop-tail configuration is Zi- t = pr? (n-i)^ • ^ n ^ ne 
other side, the partition function of a tail conformation with no loops whatsoever (i.e., of a 
nonadsorbed tethered chain) is Z t = [1% iV 71-1 . Thus the probability P^ op (l) to find a loop 
of length I next to a tail of length N — I can be estimated as P^ op (l) = oc jr+^n^zm > 
which is valid at e < e c . In the vicinity of the CAP, e ~ e c , the distribution will be given 
by an interpolation between the expressions above. Hence, the overall loop distribution 
becomes 
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The same reasonings for a tail leads to the distribution 
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In eqs.flH]) - (EJ) At, A 2 , B\, B 2 are constants. Close to CAP these distributions are expected 
to attain a U - shaped form (with two maxima at / = 1 and / ~ N), as predicted for a 



Gaussian chain by Gorbunov et al 
function for loops yields L 



131 ] . For the average loop length L the GCE-partition 

zdU(z)/dz\ z=z * = $(a - 1, [i 3 z*)/$(a, H3Z*). Close to the 



CAP, L diverges as L oc l/(e — 6c) 1 ' The average tail length S 1 can be obtain as 5 = 
Z(9Q(-2)/<92;| z=z * = — 1, /^3Z*)/[l + $(/3, ^32*)]. Again, using the properties of the polylog 
function, one can show that close to e c the average tail length diverges as S oc l/(e — Cc) 1 ^. 
Note that this behavior corresponds to a length of adsorption blob g oc l/(e — fc) 1 ^. 

i?ofc 0/ interacting loops and tails - Consider the number of configurations of a tethered 
chain in the vicinity of the CAP as an array of loops which end up with a tail. Using the 
approach of Kafri et al. [3] along with Duplantier's graph theory of polymer networks, 
one may write the partition function Z for a chain with M building blocks: M — 1 loops and 
a tail. Consider a loop of length M while the length of the rest of the chain is K, that is, 
M + K = N. In the limit of M > 1, K > 1 (but with M/A" < 1) one can show Q that 
Z ~ /if M 7 ^" 7 ^-! /if A"^-!" 1 where the surface exponent 7^ = 2 - A/> + 1) + a x + o\ 



and oi, a{ are critical bulk and surface exponents [9J. The last result indicates that the 
effective loop exponent a becomes 



« = Tat-i - Tat = " + 1 ( 7 ) 



Thus, = a — l = z/ = 0.588, in agreement with earlier Monte Carlo findings [15j. One 
should emphasize, however, that the foregoing derivation is Mean-Field-like (Z appears as 
a product of loop- and rest-of-the-chain contributions) which overestimates the interactions 
and increases significantly the value of a, serving as an upper bound. The value of a, 
therefore, is found to satisfy the inequality 1 — 711 < a < 1 + u, i.e., depending on loop 
interactions, 0.39 < <p < 0.59. 

Adsorption under detaching force - Using the GCE approach now we treat the case 
of self-avoiding polymer chain adsorption in the presence of pulling force, thus extend- 
ing the consideration of Gaussian chains by Gorbunov et al. [161 ] . Under external de- 
taching force /, the tail GCE-partition function Q(z) in eq. (pQ) has to be replaced by 
Q(z) = 1 + J2^=A(t i 3 z ) n I ' n ^\ I d 3 rP n (r) exp(,f r 1 /T) where -P n ( r ) is the end-to-end distance 
distribution function for a self- avoiding chain [17] . After some straightforward calculations 
the tail GCE-partition function can be written as 

Q(z) = l + ai ~f zfr exp(a 2 / 1 ^)) (8) 

Here the dimensionless force / = fa/ksT, the exponents ip = I — v, and 9 = (2 + t — 
35/2)/(5 - 1) with t = (/3 - 3/2 + 3i/)/(l - v) and 5 = 1/(1 - v). The function Q(z) has 
a branch point at z# = /ig 1 exp(— a2f 1 ^ u ), i.e., Q(z) ~ l/(-2* — z) l ~^ . One may, therefore, 
conclude that the total GCE-partition function has two singularities on the real axis: 
the pole z*, and the branch point It is known (see, e.g., Sec. 2.4.3. in 111 ]) that for 
iV >> 1 the main contributions to S^r come from the pole and the branch singular points, 
i.e., 



aN ~ d (z*)"^ + f( ^ y iV"^ (,#)^ (9) 

Thus, for large N only the smallest of these points matters. On the other hand, z* depends 
on the dimensionless adsorption energy e only (or, on w = exp(e)) whereas is controlled 
by the external force /. Therefore, in terms of the two control parameters, e and /, the 
equation z* (e) = z# (/) determines the critical line of transition between the adsorbed phase 



and the force-induced desorbed phase . In the following this line will be called detachment 
line. Below it, / < fo, or above, / > /d, either z* or z#, respectively, contribute to Sjy. 
The controll parameters, and /©, which satisfy this equation, denote detachment energy 
and detachment force, respectively. On the detachment line the system undergoes a first- 
order phase transition. The detachment line itself terminates for fu — > in the CAP, e c , 
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FIG. 1: (a) The 'order parameter', n, against the surface potential, e, for various pulling forces. 
The chain has length iV=128. (b) Variation of n with the pulling force, /, for several surface 
potentials. 



force fz) vanishes as fo ~ (e — e c ) u ^. This first order adsorption-desorption phase transition 
under pulling force has a clear dichotomic nature (i.e., it follows an "either - or" scenario): 
in the thermodynamic limit iV — ► oo there is no phase coexistence] The configurations are 
divided into adsorbed and detached (or stretched) dichotomic classes. Metastable states are 
completely absent. Moreover, the mean loop length L remains finite upon detachment line 
crossing. The average tail length S, on the contrary, diverges close to the detachment line. 
Indeed, at / < fo the average tail length is given by S = f 6 <&{jp — 1, z*(w)/z # (f))/[l + 
ai$(?/', z*{w) / z*{f))}. At the detachment line, z* = z#, it di verg es as S oc /d/(/d — /)• 

Reentrant phase behavior - Recently, it has been realized [l8| that the detachment line, 
when represented in terms of dimensional variables, force fr, versus temperature T, goes 
(at a relatively low temperature) through a maximum, that is, the desorption transition 
shows reentrant behavior! Below we demonstrate that this result follows directly from 
our theory. It can be seen that the solution of eq.(T5]) at large values of e (or, at low 
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FIG. 2: Plot of the critical detachment force fjj = fa/ksT against the surface potential e/ksT. 
Full and empty symbols denote MC and theoretical results. A double logarithmic plot of fo against 
e — e c with e c = 1.67 is shown in the inset, yielding a slope of 0.97 ± 0.02, in agreement with the 
prediction fjj oc (e — e c ) u ^. Shaded is shown the same phase diagram, derived by numeric solution 
of z* = z#, which in dimensional / (right axis) against T (top axis) units appears reentrant. 

temperature) can be written as z* pa e~ e /^ so that the detachment line, z* = z#, in 
terms of dimensionless parameters is monotonous, fo oc [en — hi (^3/^2 Note, however, 
that the same detachment line, if represented in terms of the dimensional control parameters, 
force fjj versus temperature Tjj (with a fixed dimensional energy e ) > shows a nonmonotonic 
behavior fp = Td[eo/Td — In(/i3//i2)] 1/ /a- This curve has a maximum at a temperature 
given by T% ax = (1 - u)e / ln(// 3 //x 2 ). 

Monte Carlo Simulation - We have investigated the force induced desorption of a polymer 
by means of extensive Monte Carlo simulations using a coarse grained off-lattice bead-spring 
model [3| of a polymer. Fig. [Th shows the variation of the order parameter n (average 
fraction of adsorbed monomers) with changing adhesive potential e at fixed pulling force 
whereas FigJTjD depicts n vs. force fa/T for various e. The abrupt change of the order 
parameter is in close agreement with our theoretical prediction. Using the threshold values 
of fi> and for critical adsorption/desorption in the thermodynamic limit iV — *• 00, one can 
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construct the adsorption - desorption phase diagram for a polymer chain under pulling shown 
in Figj2] which is among the central results of this work. The detachment lines, obtained 
from MC data and the numerical solution of z* = z& almost coincide, and the slope of fjj 
vs (e — e c ) is close to unity, according to the prediction f D oc (e — e c ) u ^. Also indicated by 
the shaded area in Figf2] is the reentrant image of the same phase diagram, obtained when 
the numerical solution of z* = 2* is plotted in dimensional units / versus T. In Fig. [3^ 




S io' 



(b) 



fa/k D T=0.0 



♦ ♦ ♦ 



, > # \ 




% s 

♦ i s 



0.01 0.02 0.03 0.04 
,1/N 



• N=128 
■ N=64 

♦ N=32 



(8-8)/k R T 



FIG. 3: (a) Tail length distribution P(s) for different surface potentials close to e c in a polymer of 
length ./V = 128 with no pulling force. In the inset P(s) at e = e c (symbols) is compared to the 
prediction, Eq. ^ (full line), (b) Average tail length S against (e — e c )/kBT plotted for various 
chain lengths in log- log coordinates. The slopes of these curves are plotted against 1/N in the inset 
and extrapolate to 1/0 in the thermodynamic limit N — > 00. 



we show the PDF of tail length at different strength of adsorption in the absence of pulling. 
This confirms the U - shape of P(s) predicted by eq.([6]). While for s — > 1 the agreement 
with eq. [6] is perfect, for s — > N long tails are slightly overestimated by eq.Q. This small 
discrepancy reflects the dominance of our "single loop & tail" approximation - multiple loops 
would effectively reduce the tail size. Fig. [3b shows the divergency of S close to the critical 
point e c . For chain of finite length N, the tail length divergence at e —>■ e c is replaced by a 
rounding into a plateau since S — >• iV but away from e c the measured slope extrapolates to 
the theoretical prediction S oc l/(e — t c ) 1 ^- In the presence of pulling force one observes a 
remarkable feature of the order parameter probability distribution - an absence of two peaks 
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FIG. 4: Distribution of the order parameter n for a pulling force /a/fcgT = 6.0 an different 
strengths of adhesion e/ksT. The chain length is N = 128 and the threshold value of the surface 
potential for this force is ed ~ 6.095 ± 0.03. The values e/A^T = 6.09 and e/ksT = 6.10 are on 
both sides of the detachment line, cf. Fig. [2j 

in the vicinity of the critical strength of adsorption, « 6.095 ± 0.03, which still keeps the 
polymer adsorbed at pulling force fa/ksT = 6.0 - Fig. HI At eo the distribution H{n) is 
flat, indicating huge fluctuations so that any value of n is equally probable. Close to eo, 
one observes a clear maximum in the distribution H(n), indicating a desorbed chain with 
n w 0.01 for e = 6.05, or a completely adsorbed chain with n « 0.99 for e = 6.15. This 
lack of bimodality in the H(n) manifests the dichotomic nature of the desorption transition 
which rules out phase coexistence. 

In conclusion, we have demonstrated that a full description of the force induced polymer 
chain desorption transition can be derived by means of the GCE approach, yielding the 
average size and probability distribution functions of all basic structural units as well as their 
variation with changing force or strength of adhesion. The detachment transition is proved 
to be of first order albeit dichotomic in nature thus ruling out phase coexistence. The critical 
line of desorption, while monotonous when plotted in dimensionless units of detachment force 
against surface potential, becomes non-monotonous in units of force against temperature, 
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thus outlining a reentrant phase diagram. In addition, we show that the crossover exponent, 
0, governing polymer behavior at criticality, depends essentially on interactions between 
different loops so that 0.39 < < 0.59. All these predictions appear in very good agreement 
with our Monte Carlo computer simulation results. 
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